Effects of Disorder on Ferromagnetism in Diluted Magnetic Semiconductors 
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We present results of a numerical mean field treatment of interacting spins and carriers in doped 
diluted magnetic semiconductors, which takes into account the positional disorder present in these 
alloy systems. Within our mean-field approximation, disorder enhances the ferromagnetic transition 
temperature for metallic densities not too far from the metal-insulator transition. Concurrently, the 
ferromagnetic phase is found to have very unusual temperature dependence of the magnetization as 
well as specific heat as a result of disorder. Unusual spin and charge transport is implied. 



Following the discovery of a ferromagnetic transition in 
Gai_a;Mn2;As at temperatures in excess of 100 K 
well above those found in counterparts based on II- VI 
semiconductors Q], there has been a surge in interest 
in the magnetic properties of diluted magnetic semicon- 
ductors (DMS). Theoretical models abound to explain 
the ferromagnetism [^-0. While it is generally accepted 
that the ferromagnetism is due to an effective interaction 
between the magnetic ions (Mn) mediated by mobile car- 
riers (holes, since Mn, a group II element substitutes for 
Ga, a group III element) , different models differ in detail, 
e.g. whether the interaction is RKKY or not, and also 
the approximations used to model the system. 

In nonmagnetic doped semiconductors, such as phos- 
phorus doped silicon there has been no evidence for 
ferromagnetism due to carriers. Indeed, carrier hopping 
at low doping concentrations in the insulating phase is 
known to induce antiferromagnetic interactions between 
localized states, leading to a valence-bond-glass like state 
down to the lowest observable temperatures In con- 
trast, ferromagnetic tendencies were detected in doped 
diluted II- VI magnetic semiconductors already in the 
insulating regime at low temperatures [Q, and subse- 
quently ferromagnetism was observed in both II- VI and 
III-V semiconductors at metallic doping densities. 

In insulating DMS, the presence of Mn has been shown 
[0,^ to overwhelm the antiferromagnetic interaction 
between charge carriers, leading to an essentially ferro- 
magnetic ground state. Monte Carlo simulations jl^ for 
II- VI DMS in the insulating phase show that the ferro- 
magnetic phase is very unusual, with a highly inhomoge- 
neous magnetic profile, leading to unconventional prop- 
erties such as M(T) curve that is not described by expan- 
sions around the critical point (critical point theories) or 
zero temperature (spin wave theories) over most of the 
ferromagnetic phase. By contrast, theoretical models for 
the metallic regime have been based on the homo- 
geneous electron gas, with a few exceptions, such as the 
possibility of phase separation [Q. 

It is well-known in conventional doped semiconductors 
that the carrier wavefunctions are those derived from an 
impurity band, for densities in the vicinity of the metal- 



insulator transition (MIT), up to factor of 3-5 above the 
MIT density, ric . Density Functional calculations jl^] 
for a lattice of hydrogen atoms show this clearly. The 
variation of critical density with uniaxial stress in doped 
Si and Ge jl^ is in agreement with calculations based 
on impurity band wavefunctions jl^. Local moments 
are known to dominate the low temperature behavior in 
doped semiconductors well into the metallic phase [p^,p |, 
and the effect is enhanced for compensated systems [[19 1. 
Raman measurements [ pO| in doped Si and infrared spec- 
tra in GaAs also show features in the metallic phase 
characteristic of the impurity wavefunction. 

Ferromagnetism is found in Gai_2,Mna;As not far from 
a MIT, with insulating behavior seen both at low and 
high Mn concentration p2]. Further, the system is 
heavily compensated p2-G4] with a carrier density only 
around 10% the Mn density. (Mn is nominally an accep- 
tor in GaAs and is expected to donate one hole per Mn). 
The vicinity of the MIT, and the large compensation, 
which implies large disorder, motivate studying a model 
that takes into account disorder as well as the impurity 
potential at the outset, to see what their effects are on 
the magnetic properties of the system. 

Given the added complications of disorder, we study a 
model based on an impurity band of hydrogenic centers 
with spin-1/2 (instead of the more complex s=3/2 wave- 
functions appropriate for acceptors), coupled to localized 
Mn d-electrons in their S=5/2 ground state. The s=3/2 
case, while technically more complicated, should yield 
qualitatively similar conclusions. The impurity band is 
described in terms of a tight binding Hamiltonian of 
the ground state impurity wavefunctions at the impurity 
sites, which are distributed randomly on the Ga sites . 
As in previous work [ pl]JT^ , the carriers are coupled to 
the Mn spins by an antiferromagnetic (AFM) Heisenberg 
exchange interaction. The Hamiltonian we study is thus: 
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Here, Ri {i = l,Nd) denotes the random positions of 
the Mn impurities, and c}^ is the creation operator of 
a hole with spin a in the bound state associated with 
the ith Mn impurity. The first term in Eq. (|l|) de- 
scribes the hopping of holes between various site, with 
the hopping matrix tij = t{\Ri — Rj\) given by t{r) ~ 
2 (1 + r/fls) exp (— r/aB)Ry 0, where the Ry is the 
binding energy of the hole, Eb, and as = eh'^/m*e^ is 
the hydrogenic Bohr radius. The second term is the AFM 
interaction between Mn spins S(i) and hole spins. Since 
Mn spins are strongly localized, the exchange integral is 
simply given by = Jexp {—2\Ri — Rj\/aB), reflecting 
the probability of finding the hole in the impurity state 
around j on the zth Mn spin. The last line in Eq. (|l|) 
describes interactions with an external magnetic field H. 

We study finite size lattices containing simple cubic 
unit cells (lattice constant a) of the zinc-blende structure. 
Nd of the Ga FCC sublattice are substituted at random 
by Mn, leading to a Mn concentration umu = 4a;/a^, 
where x = Nd/AL^. The total number of holes is Nh = 
pNd, implying a hole concentration Uh — priMn- In all 
simulations presented in this paper we choose L such that 
for the corresponding x and p, we have Nh > 50 and 
Nd > 500, so as to minimize finite size effects. Thus, 
in the absence of external magnetic fields, the problem 
can be scaled in terms of four dimensionless parameters: 
J/Eb, as /a, riha^ and x. 

In this paper, we use parameters believed to be appro- 
priate for Gai_a;Mna;As: lattice constant a = 5.65A, hole 
binding energy Ef, — 112 A meV=lRy ||2^, with a conse- 
quent Bohr radius (in our model) of = 7. 8 A |2^, and 
an exchange integral J = 15 meV p^ ]. Typical values 
of the Mn and hole concentrations are x — 0.01 — 0.05 
and p = 5 — 10% [^,^. A more comprehensive study, 
including a number of effects left out of this model, is be- 
ing completed Q. With these parameters, the typical 
hopping parameter is t(4as) = 20meV, though it should 
be emphasized that Uj are distributed over a wide range 

We treat the AFM interaction within the mean-field 
approximation (MFA), which leads to the replacement 
S{i)sj {S{i))sj + S{i){sj} - (5(«)}(sj), where = 
Cj^-^cTapCjfs. The charge carrier Hamiltonian now con- 
tains the hopping term and an effective on-site interac- 
tion Ejea/3(j)c]aCj73, with e„/3(j) = ^J2^J^j^cp{S{i)), 

and can be numerically diagonalized for any configura- 
tion {S{i)) of the Mn spins, allowing the calculation of 
the charge carrier spin expectation values (sj). In turn, 
these allow us to compute the new expectation values 
for the Mn spins, \{S{i))\ = BsiPH,),S{i) \\ H,, where 
Hi = — 'Y^^ Jij{sj) and Bs{x) is the Brillouin function. 

The process is repeated until self-consistency is reached 
at each site ||2^. As usual in MFA, the symmetry to spin 
rotations is spontaneously broken and the expectation 
values are non-zero for some direction, which we choose 
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FIG. 1. The average Mn spin SMn(> 0) and the average 
spin per hole Sh{< 0) for a typical random Mn distribution 
(full lines) and a simple cubic ordered Mn distribution (dashed 
lines), both for x = 0.00926 and p = 10% (see text). Inset 
shows the corresponding specific heats per Mn spin. 

as the z-axis (this is equivalent to having an infinitesi- 
mally small magnetic field). The average contributions 
to the total magnetization of the Mn and hole spins 
are then proportional to Smu = I ^ dYln^S^ ij)) and 
s/i = IjNh'Y^^l^sl). As might be expected, below a tem- 
perature Tc, the system develops non zero expectation 
values for the Mn and hole magnetizations, through hole- 
induced alignment of the Mn spins. In Fig.|l| we show the 
average Mn and hole spins as a function of temperature, 
T, for a system with x = 0.00926 and p = 10%, for a 
typical random Mn distribution (full lines) . For compari- 
son, we also show the corresponding results (dashed lines) 
for a system with the same Mn concentration, but with 
Mn ions arranged on a simple cubic lattice, with a (su- 
per)lattice constant oz, = a/(4a;)^/'^ = 3a. Due to their 
AFM interaction, the two expectation values have oppo- 
site signs, with the Mn spin saturating at | and the hole 
spin saturating at — i at low temperatures. (The total 
magnetization M{T) of the system has a T-dependence 
similar to that of Smu, since Mn spins outnumber holes 
ten to one, and also have a higher moment). 

The first observation is that the magnetization of the 
disordered system does not have the Brillouin-function 
shape typical for uniform ferromagnets. This is in part 
due to the small carrier density relative to the Mn spin 
density; however, an even greater effect comes from the 
wide distribution of exchange couplings and hopping in- 
tegrals, because of which many Mn spins do not order 
down to extremely low T. This is made clear in the in- 
set for Fig.|l|, which compares the specific heat for the 
two configurations on a logarithmic scale for T: the dis- 
ordered case shows a pronounced peak at temperatures 
well below Tc, and significantly below its lattice counter- 
part. Further, we find that the average Mn moment at 
T = 0.2Tc is well below the saturation value of 5/2, in 
accord with experimental data p3| , but in contrast with 
results obtained from homogeneous electron gas models. 
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FIG. 2. The average Mn spin Smk and average spin per 
hole Sh for doping concentration x — 0.00926 and p = 10%. 
In increasing order of Tc, the curves correspond to an or- 
dered, weakly disordered, moderately disordered and com- 
pletely random distributions of Mn (see text). 

A second, surprising result is that randomness in the 
Mn positions leads to a significant increase in Tc- This is 
because, in the disordered system, holes prefer regions of 
higher local concentration of Mn, where they lower their 
total (magnetic and kinetic) energy, by polarizing the Mn 
spins and hopping among several nearby Mn sites. As a 
result, these regions of higher Mn concentration become 
spin-polarized at higher temperatures than in the uni- 
form system, and the resultant Tc is increased. This is 
similar to a percolation-like situation. We caution that 
this increase may be significantly overestimated in a MFA 
such as ours, since spin fluctuations between weakly cou- 
pled polarized clusters are not treated accurately. 

On the other hand, in a disordered system, the lower 
density Mn regions have a lower than average probabil- 
ity to be visited by the holes, and as a result the Mn 
spins in these regions only align ferromagnetically at ex- 
tremely low temperatures (see Fig.|l|); such an effect is 
probably well captured by our scheme. We would like 
to emphasize the fact that the holes at the Fermi energy 
at the densities studied are either itinerant, or close to 
being so. An analysis of hole wave-functions ||2^ shows 
this delocalization, along with the higher weight of holes 
in the regions of high Mn concentration. This delocal- 
ization is responsible for alignment of the polarization at 
these high temperatures (relative to the insulating sys- 
tem ll^). Holes traveling between various high-density 
regions force the alignment of Mn spins in each region to 
be the same, in order to minimize their kinetic energy. 

This picture can be checked by "tuning" the amount 
of disorder (randomness) in the Mn positions. In Fig.|^ 
we show the average hole and Mn spins curves, for four 
types of Mn distributions, with x = 0.00926 and p = 
10%. In order of increasing Tc, they are: (a) an ordered 
Mn cubic superlattice; (b) weak disorder, corresponding 
to randomly displacing each Mn in (a) to one of the 12 
nearest neighbor sites of the underlying FCC sublattice; 
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FIG. 3. The average Mn spin Smu and the average spin 
per hole Sh for doping concentration x — 0.05 and p = 10% 
(thick lines) and 30% (thin lines) for typical random Mn dis- 
tributions (full lines) and simple cubic ordered Mn distribu- 
tions (dashed lines). 

(c) moderate disorder, corresponding to a random distri- 
bution of Mn on the FCC sublattice, subject to a con- 
straint that all Mn-Mn distances are greater than 2a; (d) 
completely random distribution of Mn on the FCC sub- 
lattice. With increasing randomness, Tc increases, while 
saturation of M (T) is simultaneously pushed to lower T. 
The sensitivity to disorder suggests that carrier density is 
not the only parameter characterizing the ferromagnetic 
behavior in DMS; indeed, since original submission of 
this manuscript, changes in M(T) curves with annealing 
time have been seen experimentally |30|] . 

We find a qualitatively similar picture holds for higher 
Mn concentrations as well as higher hole densities, 
though the effects are quantitatively less. (Our model, 
which does not include band states, is likely to be less 
accurate at high densities). In Fig. ^ we show the Mn 
and hole spins for both the simple cubic superlattice and 
the random Mn distribution on the FCC sublattice for 
X = 0.05 for two different p. While Tc is again larger in 
the random system in MFA, the percentage increase is 
smaller than in the x = 0.00926 case. Increasing the hole 
concentration from p — 10% to p = 30% makes the curves 
more Brillouin-like. This is because the fluctuations in 
the local doping are smaller at higher Mn concentrations, 
and increased hole doping further reduces the width of 
the exchange distribution. 

Our model, being based on the low doping limit, likely 
overestimates the role of disorder; however, because fer- 
romagnetism in DMS is seen at low doping densities, not 
too far from the metal-insulator transition, our work does 
strongly suggest that models based on the homogeneous 
electron gas (whether mean-field , or with perturbative 
RKKY exchange g) will not correctly capture the nature 
of ferromagnetism. In particular, it casts doubt on their 
quantitative fits to the observed Tc in Gai_2;Mna;As. Our 
calculation, while including the random positions of the 
Mn dopants, leaves out disorder effects due to compen- 



sation, as well as fluctuation effects left out in the MFA. 
At higher Mn concentrations, direct Mn-Mn interactions 
(which are needed to account for spin-glass like behavior 
seen in many II- VI DMS above 20% Mn) become impor- 
tant, while for higher hole concentrations, one may have 
to include the band states in addition to the impurity 
band, and possibly Coulomb interactions between car- 
riers as well. While these will have quantitative effects 
on the results |3|], the unusual shape of the magneti- 
zation curve and thermodynamic properties is likely to 
remain at low doping, judging from the results of nu- 
merical Monte Carlo simulations for the insulating phase 
of doped DMS Local experimental probes such as 

ESR and NMR would be especially valuable in ascertain- 
ing any inhomogeneities in the magnetization and carrier 
density profile, and help uncover the nature of ferromag- 
netism in doped DMS at these low carrier densities. 

In conclusion, the nature of ferromagnetism in doped 
DMS for low doping, not too far above the metal- 
insulator transition density, is strongly affected by disor- 
der, which may, surprisingly, aid higher Tc in this regime. 
Further, by appropriate tuning of various parameters, 
one may tailor the magnetic behavior M(H,T) in a man- 
ner not possible in simple uniform magnets ||3|]. This 
versatility makes DMS ferromagnetism near the MIT a 
very interesting problem from a theoretical point of view. 
Adding to the richness are possible effects of direct Mn- 
Mn interactions in concentrated systems (which lead to 
spin-glass behavior in undoped II- VI DMS |^^), the ex- 
istence of a ferromagnetic metal-insulator transition (un- 
like conventional doped semiconductors and amorphous 
alloys), and the likely unusual electron and spin transport 
characteristics because of disorder. 
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